LPTENS-00/20 
IEM-FT-204/00 
IFT-UAM/CSIC-00-31 
FERMILAB-PUB-00 /240-T 
ANL-HEP-PR-00-101 
EFI-2000-030 

Supersymmetric CP-violating Currents and 
Electroweak Baryogenesis 



M. Carena", J.M. Moreno^ M. Quiros"^ 

M. Seco^ t and C.E.M. Wagner^'^ 

^ Fermi National Accelerator Laboratory, P.O. Box 500, Batavia, IL 60510, USA 
^Instituto de Estructura de la Materia (CSIC), Serrano 123, E-28006 Madrid, Spain 
^HEP Division, Argonne National Laboratory, 9700 Cass Ave., Argonne, IL 60439, USA 
^Enrico Fermi Institute, Univ. of Chicago, 5640 Ellis Ave., Chicago, IL 60637, USA 



Abstract 

In this work we compute the CP-violating currents of the right-handed stops 
and Higgsinos, induced by the presence of non-trivial vacuum expectation values 
of the Higgs fields within the context of the minimal supersymmetric extension of 
the Standard Model (MSSM) with explicit CP-violating phases. Using the Keldysh 
formalism, we perform the computation of the currents at finite temperature, in 
an expansion of derivatives of the Higgs fields. Contrary to previous works, we 
implement a resummation of the Higgs mass insertion effects to all orders in pertur- 
bation theory. While the components of the right-handed stop current j~ become 
proportional to the difference H2d^Hi — Hid^H2 (suppressed by A/3), the Higgsino 

currents, j~ , present contributions proportional to both H2d^Hi it Hid^H2. For 

Hi 

large values of the charged Higgs mass and moderate values of tan [5 the contribution 
to the source proportional to H2d^Hi + Hid^H2 in the diffusion equations become 
sizeable, although it is suppressed by the Higgsino number violating interaction rate 

— 1/2 

. For small values of the wall velocity, 0.04 ^ ^ 0.1, the total contribution 
leads to acceptable values of the baryon asymmetry for values of the CP-violating 
phases ^cp in the range 0.04 ^ \ smipcp\ ~ 1. Finally, we comment on the relevance 
of the latest results of Higgs searches at LEP2 for the mechanism of electroweak 
baryogenesis within the MSSM. 



'^Present address: Department of Physics, University of Virginia, 382 McCormick 
Road, P.O. Box 400714, Charlottesville, VA 22904-4714. 



1 Introduction 



The origin of the baryon asymmetry of the Universe is one of the most important open 
questions in cosmology and particle physics. It has been long understood that, in order to 
generate the observed baryon asymmetry, three requirements need to be fulfilled: the 
non-conservation of baryon number, CP-violation and the existence of non-equilibrium 
processes j2]. Interestingly enough, at temperatures above the electroweak phase tran- 
sition temperature, T^, the Standard Model fulfills these requirements. Baryon number 
violation is induced by anomalous 13] sphaleron processes Hj, which are suppressed at zero 
temperature, but whose rate grows linearly with the temperature above Tc The non- 
conservation of CP is an essential property of the Standard Model, and non-equilibrium 
processes may be obtained through the expansion of bubbles of true vacuum, which occurs 
after the electroweak phase transition. 

In spite of fulfilling all the desired properties, the rate of the CP-violating processes 
in the Standard Model (SM) is too small to induce the required baryon asymmetry jEllZI- 
Moreover, the preservation of the generated baryon asymmetry after the electroweak 
phase transition requires a strongly first order phase transition 0, with v{Tc)/Tc > 1, 
where f (T^) is the Higgs vacuum expectation value at the critical temperature T^. For 
the experimentally allowed values of the Higgs mass, this requirement is not fulfilled in 
the Standard Model 

Supersymmetric particles lead to new radiative corrections to the Higgs effective po- 
tential at finite temperature JDj-f^- Light boson fields with relevant couplings to the 
Higgs field may induce a stronger first order electroweak phase transition The 
supersymmetric partners of the top quark are the only new bosons which couple in a 
relevant way to the Higgs boson which acquire vacuum expectation value and hence play 
a relevant role in defining the strength of the phase transition H. For sufficiently small 
values of the stop masses the strength of the phase transition is enhanced 1^ . In 
order to get values of v{Tc)/Tc > 1, however, the right handed stop soft supersymmetry 
breaking squared mass parameter, mfj, should be small or even slightly negative and the 
stop mixing mass parameter = \At — tan/5| must be smaller than ~ 0.6 mg, with 
niQ the left-handed stop supersymmetry breaking mass parameter. Under these condi- 
tions, and for niq <l-3 TeV, a strongly first order phase transition may be obtained up 
to values of the lightest CP-even Higgs boson mass as high as ~ 110-115 GeV pT| I28j. 

Moreover, supersymmetric particles lead to new, relevant CP- violating sources for the 
generation of the baryon asymmetry [29 . Several computations have been performed pUj - 
|12] in recent years, showing that if the CP-violating phases associated with the chargino 
mass parameters are not too small, these sources may lead to acceptable values of the 
baryon asymmetry. In this work, we shall perform a computation of these new CP- 
violating sources in an expansion in derivatives of the Higgs background fields. Similarly 
to Ref. j|31|, we shall use the Keldysh formalism [IHl for the computation of the CP- 

^An alternative dynamics for preserving the generated baryon asymmetry has been explored in Ref. jS]. 

^ Although bottom and tau Yukawa couplings become large for large values of tan /?, the bottom 
and tau superfield couplings to the Higgs boson combination which acquires vacuum expectation value, 
$ = Hi cos (3 + H2 sin remains small, apart from an enhancement of the <i>-trilinear coupling to left and 
right sbottoms and staus, which increases the corresponding mixings, but does not lead to an enhancement 
of the phase transition strength. 



2 



violating currents at finite temperature. We improve the computation of Ref. in 
two main aspects. On the one hand, instead of computing the temporal component of 
the current in the lowest order of Higgs background insertions, we compute all current 
components by performing a resummation of the Higgs mass insertion contributions to 
all order in perturbation theory. The resummation is essential since it leads to a proper 
regularization of the resonant contribution to the temporal component of the current 
found in Ref. [HI] and leads to contributions which are not suppressed for large values 
of the charged Higgs mass. On the other hand, we consider, in the diffusion equations, 
the contribution of Higgsino number violating interaction rate from the Higgsino fi 
term in the lagrangian, F^, that was considered in our previous calculations in the limit 

r^/T oo. 

This article is organized as follows. In sections El and El we present the detailed deriva- 
tion of the CP-violating currents for the cases of right-handed top squarks {j^ ) and 
charginos respectively, by making use of the Keldysh formalism and resumming to 
all order in Higgs background insertions. These two sections deal with all the technical 
details of the computation, with the main results given in Eqs. ()2.17p . ()3.16p and ()3.18|) . 
In section IHwe present explicit, analytical, solutions to the diffusion equations and an ex- 
plicit expression for the baryon asymmetry in the broken phase after the phase transition 
in the MSSM. In section El we exhibit the results of a numerical analysis of our solutions. 
A discussion of present Higgs mass constraints is made in section (HI and in section [71 we 
present our conclusions and outlook. 



2 The squark sector 

Our aim in this section is to compute the Green functions for left-handed (^^(a;)) and right- 
handed {tji{x)) stop fields, describing the propagation of these scalars in the presence of 
a bubble wall. The bubble wall is assumed to be located at the space-time point z, 
where there is a non-trivial background of the MSSM Higgs fields, Hi{z), which carries 
dimensionful CP-violating couplings to the left- and right-handed stops. We shall use 
these Green functions to compute the right-handed and left-handed stop currents at the 
point z. The starting point is the lagrangian for the stop system: 




where Ai is the stop squared mass matrix which depends, through the Higgs background, 
on the space-time point. 

Clearly this is not a free lagrangian, since the mass matrix depends on the space-time 
coordinates, and we must identify the free and perturbative parts out of it. In order 
to make such a selection we will expand the mass matrix around the point z'^ = (r, t) 
(the point where we are calculating the currents in the plasma frame) up to first order in 
derivatives as, 

M{x)=M{z) + {x-zYM^{z) , (2.2) 
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where we use the notation Ai^lz) = dAi{z)/dz^ , and we can spht the initial Lagrangian 



as: 



Co{x) = \dM^)\' + \dM^)\' + ( tlix) t*^{x) ) M{z) 
Ant = (x- zY ( tl{x) t*^{x) ) M,iz) ( f^^^^ ^ . (2.3) 
Let U{z) G SU{2) be the matrix that diagonahzes A4{z). We can then rewrite Cq and 

2 

1=1 

Ant = (x- zY ( tl{x) t*^ix) ) W(^)^^(^) Wt(;.) ^ , (2.4) 

where mf{z), Xii.^) = 1; 2) are the eigenvalues and eigenvectors of M{z). Note that the 
description in terms of the mass eigenstates xA^) is useful so far the Higgs field variations 
are small for propagation lengths of the order of the inverse of the width of the stop fields, 
r~^. Under these conditions, namely L^F/f^ > 1, with and being the bubble wall 
width and velocity, respectively, an expansion in derivatives is justified [311. 

Now we can write down the two point Green function for the field {xi{^) X2{x)Y' 
matrix form: 

G^{x, y; z) = G{x, y; z) + J d'^w (w - zYG{x, w; z) U{z) M^{z) U\z) G{w, y; z) + . . . 

(2.5) 

where x and y are assumed to be close to z, the point at which the current is being 
evaluated and around which the expansion is being performed (|a; — z\, \y — z\ -C T"^), 
G{x, y; z) is the two by two diagonal free Green function of the stop mass eigenstates with 
masses mi{z), the trace over internal (a = 1,2) indices being understood in Eq. ()2.5|) . 
Explicitly, the free Green functions for each of the two stop eigenstates can be written 
as H^: 

Gf = [e{p') + fs] {pt - pn 

Gf = [0{-p') + fs] {P.^ ~ Pr) 

Gf = -Pr + fB {P^ - P-) , (2.6) 

where /b = "n-^dp"!) is the Bose-Einstein distribution function, which contains the depen- 
dence on the temperature T, 

p± = \ (0 71 

pl-pi -mi{z)±2iTj\p^\ ' ^ ' ' 

and is the stop width which can be taken to be agT independently of the stop 

mass eigenstate. 
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Since we need to calculate the CP-violating currents induced by the right-handed stop 
states, we should first go to the weak eigenstate basis. The Green functions in the weak 
eigenstate basis can be obtained from the ones given above, which were computed in the 
basis of mass eigenstates, by the following expression 

G\x,y-z) = U\z)G\x,y-z)U{z) . 

Therefore, the current for right-handed stops takes the form: 

dG\x,y-z) 



lim Tr 



d{x - y)^ 



{21 



where P2 = (o"o — cr3)/2, cij being the two by two Pauli matrices and ctq the two by two 
identity matrix, is a projection matrix which allows to separate the current induced by 
the right-handed stops from the one induced by the left-handed stops. Nevertheless, since 
baryon number is conserved at this point, the total CP-violating currents induced by left- 
and right-handed top squarks must be zero, TT[d^G^{x,y)] = 0, and therefore 



-- lim Tr 

2 x,y'^z 



dG\x,y;z) 
d{x - y)f. 



(2.9) 



After integrating over the w space-time variable, and going to momentum space, we can 
write the current in terms of free Green functions of the mass eigenstates at the point z: 



P^Tt [a^U^z) G^{p; z) U{z) M^{z)U\z) G{p; z)U{z) 



as U\z) G{p; z) U{z) M^{z) U\z) G^ip; z) U{z)\ 



(2.10) 



since the contribution induced by the linear term in z in Eq. ()2.5|) trivially vanishes 
because G{p; z) only depends on |p| and We are using the notation G"{p]z) = 
dG{p; z)/dpu- Note that in the above expression only off-diagonal terms of the derivatives 
of the mass matrix M.^{z) at the point z give a non- vanishing contribution. We shall de- 
note by M.y{z), the matrix containing only the derivative of the off-diagonal terms of the 
matrix M.{z). 

The current could be simplified a little bit more by using: 



U\z)DU{z) = aiDai + -i:i[U{z)]i:i[Da^]U\z)a:i 
where D is a diagonal matrix. Then j~ iz) can be written as: 

d'^p 



(2.11) 



j^z) = - -Tt[U{z)]Tt M.iz)U{z) 



---Tt[U{z)]Tt ^M^{z)U{z) 



{2n: 



^p^'Tr [a,G{p; z)a2G\p- z)] 



d'^p 
(2^ 



p^e'^G,{p-z) GUp;z) . 



(2.12) 
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Expanding Gi in terms of one gets: 

, . 8C 



J 



371 



-Im 



Ztt 



2 CO 



Im 



TT 



271 '-^ ' 



dp 



f'B 



mk iz 



Pt{p;z)P^{p- z) 



271 m\ {z) — 1102 

where is the derivative of fs with respect to its argument and is given by 



C, = iml iz)-mliz))Ti[U{z)]Ti 



M^{z)U{z) 



(2.13) 



(2.14) 



Using now the particular value of the squared mass matrix M. for the stop system, 

ml + hi Hl{z) ht {AH,{z) - fi*Mz)) 



M(z) 



ht (AlH^iz) - fi,H^{z)) 



ml + hi Hi (z) 



(2.15) 



where ht is the top-quark Yukawa coupling, At the left-right stop mixing parameter, and 
fic the complex Higgsino mass parameter, defined as /Xc = /^exp(2V9^), with /i real (positive 
or negative). In the above, we have neglected corrections 0{g'^). In this approximation, 
the above constant vector C^, Eq. ()2.14|) . can be written as: 



= 2hl Im(A/ic) {H2{z)H'i{z) - H,{z)H!i{z)] . 



(2.16) 



Hence, in order to compute the CP-violating currents induced by the stop fields, the 
momentum integrals should be performed. Due to the form of the free Green functions, 
Eq. ()2.6|) . the integral over the temporal component of the momentum cannot be per- 
formed by standard integration methods in the complex plane. It is therefore better to 
perform the integration over the spatial components of the momentum and express the 
results as an integral function over p^, which admits a simple physical interpretation. In 
order to perform the integrals of the spatial components of the momentum, one should 
note that all functions depend only on |pp. Therefore, the angular integration can be 
trivially performed and the integral over the modulus \p\ of a function J^(|p|) can be writ- 
ten as half the integral on the whole real plane of the function J-'{x), with T{x) = JF(— x). 
Doing this, we can perform the spatial momentum integrals in Eq. ()2.13|) by means of 
standard techniques of integration in the complex plane and the residues theorem, and 
we can cast the resulting currents as: 



7- 

where 



hi Im(A/^c) {H2{z)H^{z) - H^{z)H^{z)} {Tb{z) + 5'^°6^b(^)} (2.17) 




Z2 



mi + AiVj 



(2.18) 
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and Zi is defined as the pole of P.^, i.e. 



z,{p') = y/pO(pO + 2ir^)-m?(z) (2.19) 
with positive real and imaginary parts satisfying 'R.e{zi) — /lm.{zi). 

3 The chargino sector 

For the case of the charged gaugino-Higgsino system we will follow similar steps as the 
ones we performed before to compute the stop current. In this case the role of right- 
handed stops is played by the (left- and right-handed) Higgsinos. The starting point is 
the lagrangian: 



C{x) = h,{x)d,Yhc{x) + W,{x)d,rWc{x) + (Wcix) K{x)) M{x) ( ^''^^^ ) (3.1) 



where 



hi J \W- 

From the structure of the chargino mass matrix we can write the lagrangian in the fol- 
lowing form: 

C{x) =iljR{xya^d^'^pR{x)+tpLix)^^^d''tljLix) 

+iPr{xVM{x)Mx) + MxVm\x)iPr{x) (3.2) 

where in this expression we have used 

= [It) ' = (ir 

Expanding the masses around the point z and splitting the lagrangian into a free and 
a perturbative part, we can write, to first order in derivatives: 

j0.o{x) =^Rixyaf,d^ijRix)+Mxy^^.d^i^L{x) 
+^R{xyM{z)Mx) + Mx)^ {z)^r{x) 

Cint{x) ={,x - zr {MxyM^{z)Mx) + Mx)^Ml{z)Mx)} ■ (3.3) 

Like for the scalar case we will diagonalize M{z) by means of the matrices U{z), 
V{z) e SU{2). Additional phase redefinition can be performed in order to bring the mass 
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eigenstates to be real and positive. In general, the lagrangian can be written as: 
jO-oix) =ipR{x)'^a^,d^ipR{x) + ipL{x)^af,d^ipL{x) 

Ant(x) =(x - zr {vR{xyU{z)M,{z)V\z)^L{x) + ^ L^x^V (z) MI{z)U\z)v r{x)} (3.4) 
where mi{z) are the eigenvalues of M{z) and 

ipR{x) = U{z)^r{x), LfLix) = V{z)ipL{x) 

are the mass eigenstates at the point z. 

At this point we can write the Green functions describing the propagation of the right- 
and left-handed fermion (p fields, S^^ and S^^, respectively, as 

S;''{x,y;z) =S^\x,y;z) 

+ j d^w{w- zYS^^{x,w]z)U{z)M^{z)V\z) S^^{w,y;z) 

+ j d'^w{w- zYS^^{x,w-z)V{z)Ml{z)U\z)S^^{w,y-z) 
S^^\x,y-z) =S'^^{x,y-z) 

+ j d'^w{w- zYS^^{x,w-z)V{z)Ml{z)U\z)S^^{w,y-z) 

+ j d^w{w - zYS^^{x, w- z)U{z) M^{z)V\z) S^^{w, y; z) (3.5) 

where S*^^, S^^, S^^ and S^^ denote the left-left, right-right, left-right and right-left 
Green functions of free fermions with mass mi{z). In the approximation where both 
fermionic widths are equal, we can rewrite the free fermionic Green functions in terms of 
bosonic ones as: 



^^^(p; z) =a,p^Gip; z) S^'^p; z) = '"'^'^ _ " J G(p; z) 



mi{z) 
m2{z) 

S^\p- z) = {^*f^ J^^^^ G{p; z) S'^^p- z) =a,p^G{p; z) (3.6) 

where the free Green functions G(p; z) are given by ()2.fi|l with fB^fF = —npHp^l), np 
being the Fermi-Dirac distribution function, mi{z) \mi{z)\, and ~ a^T. 

Using the relations between Green functions in the mass and weak eigenstate basis, as we 
did in the stop case, we obtain in the weak eigenstates basis, 

S^^ip; z) =U\z)S^^{p- z)U{z) S^\p- z) =U\z)S^\p- z) V{z) 

Si^{p;z) =V\z)Si^ip;z)U{z) S^\p; z) =V\z)S^\p; z)Viz) . (3.7) 



8 



The Higgsino currents can now be defined as: 

jj (z)= lim {TT[P2a''S^''{x,y;z)]]±TT[P,a''S!^'^{x,y;z)]]} (3. 



where P2 is the projection operator used in Eq. ()2.8p . 

By replacing ()3.6|) and ()3.7|) in these currents, and taking into account that, as hap- 
pened in the stop case, the contribution of the linear term in z in Eq. ()3.5|1 is zero by 
symmetry reasons, one gets El 

fnj')=l J ^{p'^Tr[a3(Wt(^)G(p;^)W(^)M,(z)Mt(;.)Wt(;.)G^(p;^)W(z) 

- U^{z)G^{p; z) U{z)M{z)Ml{z) U^{z)G{p; z) U{z) 

±V\z)G{p- z)V{z)Ml{z)M{z)V\z)GP{j); z)V{z)) 

T V\z)GP{j)] z)V{z)M\z)Mp{z)V\z)G{p- z)V{z))] 

+ Tr [U\z)G{p; z)U{z){M^'{z)M\z) - M{z)M^\z))U\z)G{p; z)U{z) 

± V\z)G{p; z)V{z){M^'\z)M{z) - M\z)M^'{z))V\z)G{p; z)V{z)\ ] . (3.9) 

The chargino mass matrix is given by 

M{z)=( ^^^"M (3.10) 



ui{z) /i, 

where we have defined Ui{z) = gHi{z). The diagonahzing matrices are jHSj 



A + A M2 ui + /i* U2 



2A{A + A) \ -iM2Ui + fi^U2) A + A 



y _ 1 f A + A M2M2 + /icUi . . . 

^2A(A + A) V -(M2W2 + /i>i) A + A ^ 



where field redefinitions have been made in order to make the Higgs vacuum expectation 
values, as well as the weak gaugino mass M2, real, 

A ={M^ - \fxc\'^ -ul + ul)/2 
A ={M^ - l^eP -ul + ul)/2 

A = (A2 + |M2Ui + /i>2|')'^' , (3.12) 
and the mass eigenvalues are given by 

(A + A + uliz)) M2 + ui{z)u2{z)fi* 



■mi{z) 



v/(A + A)(A + A) 



iA + A-uliz)),.-u,iz)u2{z)M2 
v/(A + A)(A + A) 



^Notice that the phases (fi of the mass eigenvalues, mi{z) =^ \mi{z)\ exp{iipi{z)} can be absorbed in a 
redefinition of the matrix V(z), as V(z) diag(exp{i(y9i(z)}, exp{i<y92(^)})V(z). As required, the currents 
()3.9|) do not depend on this phase redefinition. 
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Using these expressions, and the property 2A = |mi(z)p — |m2(-2)P, we can cast the 
Higgsino currents in the following general form: 

J'S, = ^"'^^'^'^^ {k(^K(^)-^i(^K(^)] I 0^,p^G,{p;z)e'^G';{p;z) (3.14) 
where e^^ = +1, and 

|(G2 + Gi)(G2 - G,)]" + f^^) (G2 - Gi) (G; - G'oj . (3.15) 

Notice that while the first term in is similar to the squark current (it is 
proportional to U2{^z^u\{^z) — U\{^z)u^2{.A\ ^^e second term in and the current 
have no counterpart in the scalar sector. The contribution proportional to U2i^z^u\{^z^ — 
U\{z)vy2{z) is proportional to the variation A/3 of the angle /? along the bubble wall. Since 
A/9 < 10^^, the corresponding contribution is suppressed. The contribution proportional 
to U2{z)u\{z) + Mi(z)-U2(-2), instead, is not affected by this suppression factor, although it 
is suppressed, for large values of tan/5, as 1/ tan/5. 

Now the integration over the spatial components of the momentum can be performed 

as in the previous section and the final currents can be cast as follows. For the current 

j'i one obtains, 
?t+ 

3t^{z) =2Im(M2,uJ {K(z)<(z) -ni(^)n^(^)] {Tf{z) ^ b^'^Qp^z)] 

+ {u\{z)-u\{z)\ K(^K(^)+«i(^K(^)]Hf(^)} (3.16) 

where the functions Tp-, Qf are defined in ()2.18j) after changing fp, fni{z) — > |mj(2;)| 

and ^ Tj^, and 

np{z) = —Re dp\l + 2fp) (^—) (3.17) 

Stt^ Jo Z2 \zi + Z2J 

with Zi{z) defined in ()2.19j) . after changing mi{z) \mi{z)\ and Fj- — > F^. Note that, 
being proportional to U2{z)—u\{z) = —u'^{z) cos 2j3{z), the second term of vanishes at 

the lowest order in the Higgs field insertions, in agreement with our previous results |i3lj . 
and it also vanishes in the case tan/5 = 1. 
For the current ?' ~ one obtains, 

j!^_{z)= 2lm{M2fic)[u2{z)ut{z)+ui{z)u^{z)]{lCF{z) + 2 [A + A] Hpiz)} (3.18) 
where the function ICp is defined as, 

)Cp{z) = --^^Re r dp' {1 + 2 M^f^-) . (3.19) 
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Note that the current ()H.18|1 appears to leading order in the Higgs mass insertion and it 
is not suppressed by A/3. However, is suppressed for large values of tan (3 (which are 
needed to push the Higgs mass beyond the most recent LEP bounds, as we will discuss 
in section ^ and, moreover, its effects on the corresponding Higgs density are damped 
by the presence of the Higgsino number violating interaction rate F^. Accordingly, its 
contribution to the BAU is small. 

A similar calculation for the neutral gaugino-Higgsino system would involve diago- 
nalization of the four-by-four neutralino mass matrix, making the analytic resummation 
treatment much more involved than for the chargino case. The analysis performed in 
Ref. jSI], to lowest order in the Higgs mass insertions, showed, as expected from a naive 
counting of degrees of freedom, that the neutralinos contribute to the Higgsino current as 
half the chargino contribution, with a total effect given by 3/2 that of the chargino. After 
resummation of the neutralino sector it would be reasonable to expect a total contribution 
to the Higgsino current equal to ~ 3/2 that of the chargino sector. However, and to be 
as conservative as possible in our calculation of the baryon asymmetry, we would just 
consider as source of baryon number the chargino current (remember that left-handed 
squarks are assumed very heavy and decouple from the thermal bath) that was computed 
in this section, keeping in mind that an enhancement factor ~ 3/2 might appear after a 
rigorous calculation of the currents in the neutralino sector. 

4 The baryon asymmetry 

To evaluate the baryon asymmetry generated in the broken phase we need to first compute 
the density of left-handed quarks and leptons, ni, in front of the bubble wall (in the 
symmetric phase). These chiral densities are the ones that induce weak sphalerons to 
produce a net baryon number. Since, in the present scenario, there is essentially no lepton 
asymmetry, the density to be computed in the symmetric phase is til = uq + XlLi ^Qi 
where the density of a chiral supermultiplet Q = (g, q) is understood as the sum of 
densities of particle components, assuming the supergauge interactions to be in thermal 
equilibrium, rig = + rig. If the system is near thermal equilibrium, particle densities, 
TT-j, are related to the local chemical potential, /Xj by the relation rij = /cj/ijT^/6, where ki 
are statistical factors equal to 2 (1) for bosons (fermions) and exponentially suppressed 
for particle masses mj much larger than T. For the calculation of the density hl we will 
use the formalism described in Refs. [SOlinil- 

We will consider those particle species that participate in fast particle number changing 
transitions, neglecting all Yukawa couplings except those corresponding to the top quark. 
In this approximation, there is no left-handed lepton number contribution to Ul. By 
introducing strong sphaleron effects, first and second family quark number is generated. 
Assuming that all quarks have nearly the same diffusion constant it turns out that |3()j . 
"f^Qi = ^Q2 = 2(^Q + nr), and then. 



riL = 5nQ + Aht . 




*We use, for the third family, the notation Q = Q3, T = T3. 
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In general we will relate particle number changing, or fermion number violating, rates 
Tx with the corresponding rates per unit volume •jx, as. 



X 



67X 
T3 



(4.2) 



The involved weak and strong sphaleron rates are: 



6 hZnnc (y.,,T, 



6 Kss - a'^T 



(4.3) 



respectively, where = 20 ± 2 [44] and = C^(l)- The particle number changing 
rates that will be considered both in the symmetric and in the broken phase are: Tyz, 
corresponding to all supersymmetric and soft breaking trilinear interactions arising from 
the htH2QT term in the superpotential, Fy^, which corresponds to the supersymmetric 
trilinear scalar interaction in the Lagrangian involving the third generation squarks and 
the Higgs Hi, and F^, which corresponds to the ficHiH2 term in the Lagrangian. There are 
also the Higgs number violating and axial top number violation processes, induced by the 
Higgs self interactions and by top quark mass effects, with rates Th and F^, respectively, 
that are only active in the broken phase. 

We will write now a set of diffusion equations involving uq, tit, (the density of 

Hi = {hi, hi)) and (the density of H2 = (^2, ^2)), and the particle number changing 
rates and CP-violating source terms discussed above. In the bubble wall frame, and 
ignoring the curvature of the bubble wall, all quantities become functions of z = r + f ^t, 
where is the bubble wall velocity. The diffusion equations are: 



v^Uq =DgnQ - F 



Y 





riT 


riH + prih 


- r 




riT 


kg 










kx 



-6F. 



2^Q riT ^^ uq + Ut 
kn kx ks 



+ 7q 



(4.4) 



v,.,nn 



--Dgu'^ + Fy 



riT 
kx 



riH + puh 



+3F, 



2^-^ + 9 



kn 

Uq + riT 



+ F, 



kn 



kn 



k 



B 



IQ 



riT 
kx 



(4.5) 



v^n'jj =Dhn'^ + Fy 



v,.,n 



--DhUl + pVy 



nQ__nT uh + puh 
kq kr 

^ _ ^ 

kg kr k 



kn 

riH + rih/p 



H 



:r. + 4F,) ^+75_ 



(4.6) 
(4.7) 



where all derivatives are with respect to z, Dq ^ 6/T and Dh ~ 110/T are the cor- 
responding diffusion constants in the quark and Higgs sectors [lEj, uh = '^^2 + ^//n 
Uh = nu-i ~ rtHi, k-H = kni + ku-^i Ty = Vy^ + Fy^ and pFy = Fyj — Fy^. The parameter 
p is in the range < p < 1. In previous analyses jSHl IHH EO] the limit F^ — 00 was 
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implicitly considered, leading to the solution rih 0- However, as we will see, for finite 
values of we obtain non- vanishing values of the density nh. 

For the sources 7q in Eqs. ()4.4|) - ()4.7|) we will follow the formalism of Refs. [301 El] 
where 7x — Jx/'^x, 'Tx being the corresponding typical thermalization time. Thus we will 
use as sources of our diffusion equations, 

^Q--v^ hi Tj \m{At^l,) H\z) I3'{z) {^b(^) + gB{z)} 
^H^-- 2v^g''T^lm{M2f^,){H\z)P'{z) [J'piz) + gp{z)] 

+ H\z) cos 2/3{z) [H{z)H'{z) sin 2(3{z) + H^{z) cos 2(3{z)/3'{z)] Hf{z) } 
~ 2v^g'^Tf^lm{M2i^c) [H{z)H'{z) sm2(3{z) + H\z) cos2(3{z)f3'{z)] 

{}CFiz) + 2{A + A)nFiz)} . (4.8) 

Notice that our sources, Eq. ()4.8p . are proportional to the wall velocity v^, and so die 
when the latter goes to zero, which is a physical requirement. 

We can find an approximate solution for riq and ht by assuming that Fy and Tss are 
fast so that ng/kQ -nr/kT - {riH + pnh)/kH = 0(1/Fy) and 2nQ/kQ - nr/kT + 9 {ng + 
i^T)/kB = C^(l/Fss). In this case we can write 



riT 



kg (Qfcr - ks) 
ku {ke + 9kQ + 9A;t) 
_ kr {9kQ + 2kB) 

kn {ks + 9kQ + 9fcr) 



{uh + puh) +0 



1 

rv 



Y 



(4.9) 



If the left-handed third generation squarks were light (mg ~ T) we could expect that all 
super symmetric and supersymmetry breaking interactions arising from the htH2QT term 
in the superpotential are in thermal equilibrium and similar in size, so that Fy^ ~ Fyj, or 
p <^ 1. In such case, which was considered in Ref. the influence of Uh in the quark 
densities uq and ut, through Eqs. ()4.9j) . is p-suppressed although this suppression can 
be arguably mild depending on the particularly chosen value of p. However, in the case 
where left-handed squarks are heavy {mg S> T), as preferred to get a good agreement 
of the MSSM with electroweak precision measurements, their corresponding interactions 
decouple, Fy^ ~ and p ~ 1. This is the case we will consider from here on. 

We now take (for p = 1) the linear combinations of Eqs. ()4.4j) . ()4.5j) . ()4.fj|l and ()4.7j) 
which are independent of Fy and F^^. They are given by. 



[uq + 2nrp - n'fj] =Dq [uq + 2 n'^] - Dh + F^ 



[u'q + 2 - n^] =D, [n'^ + 2n'i]- nl + F, 
+ [r. + 4F,]!^-(7Q + 75„ 



riQ riT 
kq kr 



(4.10) 



(4.11) 



When Uq and rix are replaced by the explicit solutions of Eqs. ()4.9p . as functions of uh 
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and Uh, Eqs. ()4.1()j) and ()4.1H) yield the system of coupled equations for uh and Uh- 



where the sources are 



with 



/± = ^ + , (4.13) 



G = kH{9kQ + 9kT + ks) , (4.14) 



and A, T> and Q are the 2x2 matrices, 



(4.15) 



with 



— _ F — _ G 

The system fj4.12|) amounts to equations for uh and n/j, with sources induced by 7g 
and 7^^, and by the same densities nH,h and their derivatives. It can be re-written as, 

u'h =Dn"H-VnH + f+ + A/+ (4.17) 

n'^ =Dh nl - + 4r^] ^ + v^Uh - Dh n'^ + + 7iy_ - (4-18) 

kh i^h 

where 

D=D^ + Dh, r = r^ + r;, (4.19) 

= - ^-^'/^ + ^« - . (4.20) 

We have solved the system ()4.12|) numerically and the results are presented in sectional 
However a very useful analytical approximation can be worked out as follows. Using 
Eq. ()4.6j) and the approximate relations ()4.9j) we can write for Uh the following equation, 

v^n', = D,nl-[T, + AT,]^ + ^^ . (4.21) 
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In this way the equation for rih has been decoupled from the other equations and can be 
easily solved. On the other hand, from ()4.17|1 and the expression for A/+, Eq. ()4.2()j) . we 
see that rih acts as a source for uh, and the equations ()4.17|) and ()4.2H) can be solved 
analytically. 

We will only quote the solutions in the symmetric phase {z < 0) since that would be 
needed to compute the baryon asymmetry from niiz), as we will see. Finally we will 
impose boundary conditions nft(±cx)) and ni^(±oo) and continuity of the functions and 
first derivatives at 2; = 0. 

From Eq. ()4.2H) we obtain the solution for nh{z), for 2; < as, 



nh{z) = Ah 



(4.22) 



where 



and 



Ah 



^vl + AT.Dh + ^v^ + AT^D, 



-C/3+ 



h JO 



,±y/vl + 4r2Dh] 



T2 



2Dh 
Th + ^T^ 



k 



H 



4r 



k 



H 



(4.23) 



(4.24) 



1/2 

Note that, from expression ()4.24j] . the coefficient Ah behaves as F^ , in the limit of large 
F^, and so the Uh density tends to zero when F^ tends to infinity, as anticipated. 
From Eq. ()4.17|) . the solution for nniz), for z < is given by 



Ah 



(4.25) 



where 



Bh =Ao + At 



F + G 



Dga+ 



D 



1 



1 



1 



- Dg{a+ + A+) + ^ 



DX 



A> 



[3- 



+ 



F + G (A+-/3+)(A+-/3_) 



[FA_ 



F(A+-/3_) 
+ D,A+) - GTm/kn] (4.26) 



and 



A 



H 



- Ah 



Dq a+ - 



F + G D 



(4.27) 
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with 



2 

1 



A==^ / dCUi^)e-'^'+ (4.28) 



and 



A± = i ± ^Jvl+ATD^^ . (4.29) 

Since we assume the sphalerons are inactive inside the bubbles, the baryon density is 
constant in the broken phase and satisfies, in the symmetric phase, an equation where ul 
acts as a source jSHI and there is an exphcit sphaleron-induced relaxation term 



v^n'j^{z) = -6{-z) [npT^snLiz) + TZnB{z)] (4.30) 
where n^? = 3 is the number of families and TZ is the relaxation coefficient 



n=-^nFT^s. (4.31) 

Eq. ()4.3()j) can be solved analytically and gives, in the broken phase 2; > 0, a constant 
baryon asymmetry, 

n^ = ^ n^r^ 1^ rf^nz.(z) e'"^'^- . (4.32) 



Using now the explicit solutions for nn and rih given in Eqs. ()4.25|) and ()4.22|) . we can 
cast the explicit solution for the baryon asymmetry as, 

ns = np —. — Ti — riTi — , n; N \ ^ 7 ( (4-33) 



ku (ke + 9kQ + Q/ct) i 7^ + v^a+ DTZ + vl 

where all symbols used in Eq. ()4.33|) have been previously defined. 

The validity of our analytical approximation is guaranteed by the dominance of nn 
over n/i, which in turn is related to the tan/5 suppression of 7^ and the presence of 
r^. In fact were we working in the limit ^ 00 we would find that the density nh is 
negligible. On the other hand, in the limit — and tan /3 ~ 1 we would really expect 
"^h > nn, due to the dominance of over 7^^, at least for large values of where the 
Ap suppression of jj^^ is more severe. However small values of tan /?, as we noticed earlier 
in this paper, are strongly disfavored in our scenario by recent LEP bounds on the Higgs 
mass. Hence, we have found that the analytical approximation is accurate with an error 
which depends on the chosen values of the supersymmetric parameters, but it is always 
much smaller than the other uncertainties involved in the final calculation. In sectional 
we will provide explicit comparison with the numerical result, while all plots will be done 
using the numerical solution of system ()4.12p . 
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5 Numerical results 



In this section we present the numerical results for the baryon asymmetry computed in 
section 13 and, in particular, of the baryon-to-entropy ratio r] = ub/s, where the entropy 
density is given by. 



45 



■9effT' 



{5.1] 



with Qeff being the effective number of relativistic degrees of freedom. The profiles H{z), 
P{z) have been accurately computed in the literature jilSlEZI- For the sake of simplicity, 
in this paper we will use a kink approximation [HI] 



H{z) =-v{T) 1 -tanh 



a 



1 - 



2z 



(3{z) =P--Af3 



1 + tanh 



a 



(5.2) 



This approximation has been checked to reproduce the exact calculation of the Higgs 
profiles within a few percent accuracy [23; provided that we borrow from the exact 
calculation the values of the thickness L^^/2a and the variation of the angle j3{z) along 
the bubble wall, A/5, as we will do. In particular we will take a ~ 3/2, L^^ = 20/T, and 
we have checked that the result varies only very slowly with those parameters, while we 
are taking the values of A/5 which are obtained from the two-loop effective potential used 
in our calculation. 

The calculation of the wall velocity is a very complicated phenomenon involving 
the hydrodynamics of the bubble interacting with the surrounding plasma. Some progress 
has been recently reported in this direction indicating that, in the case of the MSSM, 
the wall is extremely non-relativistic, v^j -C 1, and can be as slow as = 0.01. Unless 
explicitly stated, in the numerical analysis of this section, we adopt the value = 0.05, 
although the variation of the baryon asymmetry with respect to v^^ will also be analyzed. 
The widths, F^, F/i and Fy are as in Refs. [101 EI], while we are taking ~ 0.1 T 
and p = 1, in agreement with the large value we use for the left-handed third-generation 
squark masses, mq ^ 1 TeV, which makes them decoupling from the thermal bath. On the 
other hand, and consistently with the latter assumption (which is required to render the 
MSSM in agreement with the Higgs mass bounds coming from LEP), the contribution to 
ub from the squark source, 7q, is negligible. The "observable" value for rj consistent with 
Big Bang Nucleosynthesis (BBN) has been considered to be t^bbn ~ 4x 10"^^ [H]. Finally 
we will consider the third generation squark mass and mixing parameters, mg = 1.5 TeV 
and At = 0.5 TeV, and tan/? = 20 and have checked that, for all plots in this section, the 
phase transition is strong enough first order, v{Tc)/Tc ^ 1, and the Higgs mass is, within 
the accuracy of our calculations, ^ 110-115 GeV. These values are in rough agreement 
with present 95 % CL. bounds on the Higgs mass coming from LEP, or even with the 
present excess of events observed at LEP, consistent with the detection of a SM-like Higgs 
at the runs with the highest center of mass energies, ^/s > 206 GeV. We will comment 
more about the LEP constraints in the next section. 
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Figure 1: Plot of ?7/?7bbn as a function of /i for M2 = rriA = 100 GeV (thick solid 
curve), niA = 150 GeV (dashed curve) and tua = 200 GeV (dash-dotted curve), and the 
rest of parameters as indicated in the text. The thin solid curve corresponds to the case 
rriA = 100 GeV when the approximate analytical solution in ()4.33p is used. 

In Fig. [U we plot the ratio ry/r^BBN for the values of the supersymmetric parameters that 
have just been described, M2 = sin ip^ = 1 and several values of the pseudoscalar Higgs 
mass rriA- Therefore, since rj is (almost) linear in sin (f^, one can read from Fig.^the value 
of 1/ sin(^^ that would reproduce ?7bbn- This observation applies to all plots presented 
in this section, where we have fixed sin(y9^ = 1. It follows that the region of parameters 
where we find |?7/?7bbn| < 1, is forbidden in all plots. For the value itia = 100 GeV, we 
have presented both the exact result (thick solid curve), based on the numerical solution 
of Eqs. ()4.12j) . and the approximate result (thin solid curve), based on the approximate 
analytical solution ()4.33|) . We see that for values where ub/s is sizeable the discrepancy 
between the analytical and the numerical result is < 30 %. For the other curves in Fig. Q 
as well as for the rest of plots in this paper, we will use the (exact) numerical solution of 
Eqs. fB3^ . 

We are, in Fig. HI close to the resonance region discussed in Ref. [HI], which is smoothed 
by the all order resummation in Higgs mass insertions. The departure from the resonance 
is exemplified in Fig. |2l where we plot rj/rj^BN as a function of /i for M2 = 200 GeV, 
TTiA = 150 GeV and the other supersymmetric parameters as in Fig. ^ 

In Fig. El we plot ///-/^bbn as a function of tjia for M2 = /i = 200 GeV (solid curve) 
and M2 = 200 GeV, /i = 300 GeV (dashed curve), and other supersymmetric parameters 
as in Fig. [H Finally in Fig. IH we plot rj/rjBBTSs as a function of for M2 = /i = 200 
GeV, rriA = 150 GeV and the other parameters as in Fig. ^ The maximum of this curve 
comes from the interplay between the relaxation and source terms in the equation for ub, 
Eq. gM- 

The numerical results exhibited in the plots of this section are an improvement of 
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Figure 2: Plot of rj/riBB^ as a function of /x for M2 = 200 GeV and niA = 150 GeV. 

our previous results, Ref. jHI], and include the all order resummation of the Higgs mass 
insertions in the current determination, as well as inclusion of finite F^-effects in the 
diffusion equations. Since the first of these effects smooths out the resonant behaviour, 
which enhances the determination of hb for M2 = while the second one slightly 
enhances ub, our present numerical results are in rough agreement with those of Ref. \M\ . 
On the other hand if we compare our numerical results with the recent ones of Ref. |42j . 
that use WKB methods and values of p < 1 to deduce the source terms in the diffusion 
equations, we observe a discrepancy of a few orders of magnitude. However, we have 
been communicated j42j by the authors of Ref. |l2] to have detected a problem in their 
numerical codes which enhances their numerical results by some orders of magnitude and 
that might explain part of this discrepancy. As explained above, large values of mg, 
implying p = 1, are necessary in order to fulfill the present experimental Higgs mass 
bounds. 

6 Higgs mass constraints 

In this section, we shall comment on the constraints coming from Higgs searches at LEP. 
The LEP experiments at CERN have collected data during the year 1999 at various 
energies between 192 GeV and 202 GeV, for a total integrated luminosity of about 900 
pb~^. A combined limit on the Standard Model Higgs mass of about 108 GeV at the 
95 % C.L. was obtained, due to the absence of any significant Higgs signal in the LEP 
data 120] • Preliminary results of this year run j3T] show that this limit moved up by a 
few GeV (up to about 113.2 GeV). More interesting, a slight excess of events, about 3 
standard deviations above the SM predictions, has been observed, consistent with a SM 
like Higgs in the range of masses of about 113-116 GeV. 

The present Higgs mass constraints become particularly relevant for small values of 
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Figure 3: Plot of ry/r^BBN as a function of for M2 = 11 = 200 GeV (solid curve) and 
M2 = 200, /i = 300 GeV (dashed curve). 

tan/3, tan/3 < 5. In this case, due to the behaviour of the Higgs boson couplings to 
fermions and gauge bosons, the SM Higgs mass constraints translate with almost no vari- 
ations into a bound on the lightest CP-even Higgs boson massH For values of v{Tc)/Tc > 1 
and tan /3 < 5, and for left-handed stop masses smaller than ~ 3 TeV, the lightest CP-even 
Higgs mass never exceeds 105 GeV. Therefore, the mechanism of electroweak baryogenesis 
demands either values of tan (3 > 5 or unnaturally large values of mg H. 

Large values of tan (5 move the value of the Higgs boson mass, with relevant couplings 
to the gauge bosons, to larger values. However, if the values of the left-handed stop 
parameters are restricted to be below 3 TeV, for v{Tc)/Tc ^ 1, the Higgs mass cannot 
exceed 115 GeV. Observe that these values are a few GeV higher than those obtained 
previously in Ref . [2^ , since in that reference we restricted ourselves to the case of left- 
handed stop masses below 1 TeV. The observed excess of events, with bb invariant masses 
of about 114 GeV, would be consistent with electroweak baryogenesis for large values of 
tan /3 and large values of the left-handed stop mass parameters ttiq > 1 TeV, as the ones 
considered in the previous section. 

What would happen if the excess of events present at LEP would not correspond to 
a Higgs signal, but would turn out to be a statistical fluctuation with the final outcome 
of an ultimate exclusion limit for a SM-like Higgs with mass below 115 GeV? To analyze 
this, let us stress that at large values of tan (3 the coupling of this Higgs boson to bottom 

^In the presence of CP-violation, the Higgs mass eigenstates will not be CP-eigenstates. In our analysis 
we have used the CP conserving structure for the Higgs sector. This should lead to a good approximation 
if CP-violating effects in the Higgs potential are small, as happens when arg(/ic^t) — 0. A more general 
analysis, similar to the one performed in Refs. |52l I53| would be appropriate to consider more general 
CP-violating effects. 

^ We have checked that, for the values of the stop mixing parameters consistent with electroweak 
baryogenesis, no significant modification of these bounds is obtained after considering CP-violating effects 
in the Higgs potential ,52, . 
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Figure 4: Plot of ///■/^bbn as a function of Vyj for M2 = = 200 GeV and = 150 GeV. 

quarks can be significantly lower than in the SM [SH E2| with a corresponding reduction 
of the Higgs mass bound. These variations can only occur for small values of the CP-odd 
Higgs mass m^i, of order of the lightest CP-even Higgs boson mass. Unlike the case of 
small values of tan/5, for values of tan/5 > 10, the values of v(Tc)/Tc are only weakly 
dependent on the exact value of m^. Intuitively, this can be understood by the fact 
that for large values of tan/?, the CP-odd Higgs can be approximately identified with the 
imaginary part of the neutral component of the Higgs doublet Hi, while the Higgs doublet 
which acquires vacuum expectation value is mainly H2 (f2 ^ "Wi)- 

For the values of At and fi consistent with electroweak baryogenesis, a reduction of the 
coupling of the CP-even Higgs boson to the bottom quark would demand not only small 
values of ~ 100-150 GeV, but also large values of tan/? > 10 and of \fiAt\/mQ > 0.1 
(the larger tan /3, the easier suppressed values of the bottom quark coupling are obtained). 
We have checked that, assuming small CP- violating effects in the Higgs potential and in 
the Higgs-fermion couplings, and for values of mq ~ 1 TeV, v{T^/Tc ^ 1 and < 500 
GeV, a significant reduction of the coupling of the Higgs to bottom quarks only occur 
for tan j3 ^ 30. Therefore, if the excess of events observed at LEP is not associated with 
a Higgs signal, strong constraints on the electroweak baryogenesis scenario within the 
MSSM will be obtained. 

7 Conclusions 

In this article, we have performed a computation of the scalar- and fermion- CP-violating 
currents induced by the expansion of a true-vacuum bubble in the false vacuum plasma, 
within the framework of the minimal supersymmetric standard model. We made use of 
the Keldysh formalism and we have defined a systematic way of obtaining the currents 
in an expansion of derivatives of the Higgs fields, to all orders of the Higgs background 
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insertions. 

Although our method is similar to the one previously used by some of us in Ref. 
our results differ from those presented in our previous work in several respects. First 
of all, they include a resummation of corrections associated with higher order of the 
Higgs background insertions. These corrections have two important effects. The first one 
is to dilute the resonant behaviour obtained in Ref. [SI] for values of = M2. The 
second one is the appearance of a contribution proportional to H2d^Hi + Hid^H2 to the 
vector Higgsino current {z). This means that, as first observed in Refs. [1011101; the 
vector Higgsino current does not vanish for large values of the CP-odd Higgs mass. Our 
method provides a way of obtaining the value of this non-vanishing contribution in a self- 
consistent way. In addition, we have also computed the axial Higgsino current j,~ (z), 
whose components are proportional to H2d^Hi Hid^H2- Therefore, as first observed in 
Ref. [121, the chiral current is not suppressed for large values of the CP-odd Higgs mass 
and hence may become relevant in this regime. 

The vector and axial Higgsino currents, j.~ (z), were used to determine the baryon 

asymmetry of the Universe, n^/s. The computation of demands the solution of diffu- 
sion equations, with sources determined through (z). Following the method developed 

in Refs. [SHI Ell, we assumed that the sources are proportional to the temporal component 
of the currents, with a constant of proportionality given by the Higgsino width. Within 
this approximation, we computed the functional dependence of on the soft super- 
symmetry breaking parameters and on the bubble wall parameters. The most important 
parameters turn out to be the gaugino and Higgsino mass parameters, |/ic| and M2, their 
relative phase arg(/icM2), (equal to in the basis in which M2 is real) as well as the 
CP-odd Higgs mass tjia and tan (3. We have also required that the condition of preserva- 
tion of the baryon asymmetry v{Tc)/Tc 1 is fulfilled, what demands a light right-handed 
stop and, due to the present Higgs mass constraints coming from LEP (see section [HI), 
also large values of the ratio of Higgs vacuum expectation values, tan P > 5. 

Under the above conditions, we have determined the value of n^, compared to the 
value predicted by Big Bang Nucleosynthesis, for a value of simp^ = 1. The ratio of the 
theoretically obtained to the BBN predicted baryon asymmetry can be reinterpreted as 
the inverse of the value of siiup^ needed to obtain a value of ra^ in agreement with the 
BBN predictions. We conclude that, for small values of rriA — 100 GeV and |yu| ~ M2, 
values as low as 9?^ ~ 0.04 can lead to acceptable values of ub- The predicted value of the 
phase increases for larger values of rriA and/or for 7^ M2, but still there is a large 
fraction of parameter space in which the computed baryon number is in good agreement 
with BBN predictions, for phases such that siny9^ ~ 0.04-1. 

Values of ipf^ > 0.04 can lead to acceptable phenomenology if either peculiar cancella- 
tions in the squark and slepton contributions to the neutron and electron electric dipole 
moments (EDM) occur |S3], and/or if the first and second generation of squarks are 
heavy This second possibility is quite appealing and, as has been recently demon- 
strated [SIl, leads to acceptable phenomenology, including the dark matter constraints 0. 

^Third generation squarks would still contribute to the neutron and electron EDM, via two loop 
diagrams involving the would-be CP-odd Higgs boson (SSj. These contributions can become sizeable at 
large values of tan/3, although they tend to be suppressed for small values of the mixing in the stop 
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Another important observable which, similarly to the value of the baryon number, 
depends on the precise value of the mass parameters in the gaugino, Higgsino and third 
generation squark sectors, as well as on the charged Higgs mass, is the rate of the rare de- 
cay 6 — > S7 fSU' . For small values of the charged Higgs and stop masses, and for moderate 
values of At/mq and |/i|/mQ, the chargino-stop contribution, as well as the charged Higgs 
contribution, may become large for large values of tan/? ^6T H I^ H IH^. In scenarios with 
heavy first and second generation squarks, however, flavor violation couplings involving 
the third generation squarks could be non-negligible jSE] and therefore the gluino-sbottom 
contributions to this rare decay rate may be enhanced IHl] • Since these last contributions 
are strongly model dependent, and may be larger than the charged Higgs and chargino- 
stop ones, we have not imposed the 6 — > 57 constraints in our analysis. 

Finally, we have discussed the effect of the Higgs mass constraints coming from LEP. 
The preliminary data coming from the LEP experiments imply a lower bound on the 
mass of a SM-like Higgs boson of about 113 GeV. A small excess, consistent with a SM- 
like Higgs boson with a mass slightly above that value has also been observed. These 
relatively large values of the Higgs mass are consistent with electroweak baryogenesis 
within the MSSM if the value of tan/3 is large, tan/5 > 5, if the left-handed stops are 
heavy mq > 1 TeV, and if the stop mixing parameter is not small. At ^ 0.25 mq. On 
the other hand, for these values of the Higgs mass, values of > 0.4 mq make the phase 
transition weaker, leading to values of v{Tc)/Tc that are in conflict with the condition of 
preservation of the baryon asymmetry. It is important to emphasize, however, that even 
if the CP-even Higgs boson coupling to the gauge boson is SM like, it can evade the LEP 
bounds if its coupling to the bottom quark is strongly suppressed, what can occur for 
very large values of tan (3, tan f3 > 30. More relevantly, if the excess of events at LEP has 
its origins in the presence of a SM-like Higgs boson of mass of about 113-115 GeV, one of 
the predictions of electroweak baryogenesis, namely the presence of a light neutral Higgs 
boson with SM-like couplings to the gauge bosons and a mass not larger than 115 GeV 
would have been fulfilled. 
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